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NOMENCLATURE 

a  correction  foctor 

A  small  negative  number 

Aq  positive  number 

b  correction  factor 

bi  positive  n^unber  where  t  =  1, 2 
B  positive  number 

C2  heat  capacity  of  R2 

Q  constant  where  i  =  0, 1, 4 

d  unit  of  time,  day 

/  mass  flux  of  water 

F  function  defined  by  eq  A1 

I  function  defined  by  eq  12b 

/  functitm  defined  by  eq  23c,  23d,  25b  and 

25c 

Ji  function  where  i  =  1, 2, 3 

it  thermal  conductivity  of  a  frozen  fringe 

ito  fliermalconductivityoftheunfrozenpart 

of  the  soil 

ki  thermal  conductivity  of  an  ice  layer 

Kq  hydraulic  conductivity  in  flte  unfrozen 
p^  of  the  soil 

K:  empirical  function  defined  by  eq  4a 

where  i  =  1,2 

Kii  limiting  value  of  as  x  approaches  nj 
while  X  is  in  Ri,  * 

XjQ  limiting  value  of  X^  as  x  approaches  ng 

while  X  is  in  Kj,  t  =  1,2 

L  latent  heat  of  fusion  of  water,  334  Jgr^ 

m  location  of  the  free  end  of  die  column 

n  boundary  in  Rq- 

n-t  botmdary  with  i  =  0, 1  where  Hq  denotes 

the  boundary  where  T  =  0°C  and  the 
interface  between  an  ice  layer  and  a  fro¬ 
zen  fringe. 

n^o  boundary  between  Rio  and  Rii 

n*(nf)  nei^dx»hoodofniWherex<fij(x>ni) 

P  pressure  of  water 

P,  overburden  pressure 

Pq  valueofPatrto 

P„  value  of  P  at  n 

AP  defined  by  eq  10 

q  heat  flux 

Rq  unfrozen  part  of  flte  soil 

Ri  frozenfringe 

R2  ice  layer 

Rio  region  where  0  2:  r>Ti* 


Rii  regiwi where TJ*  >TiTi 

region  in  tfie  diagram  of  temperature  gra¬ 
dients  where  an  ice  layer  melts 
R,  region  in  the  diagram  of  temperature  gra¬ 

dients  where  the  steady  growfli  of  an  ice 
layer  occurs 

RS  boimdary  between  R,  and  Ru 

Ri*  boundary  between  R,„  and  R, 

R„  region  in  ttie  diagram  of  temperature  gra¬ 

dients  v^iere  die  steady  growth  of  an  ice 
layer  does  not  occur 
T  temperature 

Tq  temperature  at  kq  =  0°C. 

Ti  temperature  at  rti 

Ti  tmnperature  at  rii  at  die  formation  of  the 

final  ice  lens 

temperature  at  ni  at  the  phase  equilib¬ 
rium  of  water 

Tf  average  temperature  gradient  in  Ri 

U  defined  by  eq  29a 

Lfo  defined  by  eq  29b 

V  defined  by  eq  33a 

Vq  defined  by  eq  33b 

X  spatial  coi.  rdinate 

y  segregation  potential  function 

y**  valueofyatTi  =  Ti* 

z  defined  by  eq  24c 

Oq  absolute  value  of  the  temperature  gradi¬ 

ent  at  hq 

tti  absolute  value  of  die  limiting  tempera¬ 

ture  gradient  as  x  approaches  rii  whilexis 
in  R2,  defined  by  eq  9 

Y  constant,  1.12  MPa  ®C’^ 

S  thickness  of  a  frozen  fringe 

Sq  distance  between  n  and 

Xo  omstant  defined  by  eq  24e 

Xi  constant  defined  by  eq  24f 

o  effective  pressure  defined  by  eq  13a 

%i  empincalfuncti(mofTidefinedbyeql3b 

^1  empirical  function  of  Ti  defined  by  eq  13c 

*  superscript  used  to  indicate  the  value  of 

any  variable  evaluated  at  die  formatim  of 
die  final  ice  lens 

**  superscript  used  to  indicate  die  value  of 

any  variable  evaluated  at  tiie  phase  equi¬ 
librium  of  water 
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Dependence  of  Segregation  Potential  on  the 
Thermal  and  Hydraulic  Conditions  Predicted  by  Model  Mi 

YOSHISUKE  NAKANO 


INTRODUCTION 

We  will  consider  the  (nte-directi(»ud  steady 
growthofanicelayerimderanovetburdenpies- 
stue  P,  (2;  0).  Let  ttte  freezing  process  advance  from 
the  top  down  and  the  coordinate  x  be  positive 
upw2uxls,  with  its  origin  fixed  at  some  point  in  file 
unfixizen  part  of  file  soil.  A  freezing  soil  in  this 
problem  may  be  considered  to  consist  of  three 
parts:  the  imlFrozen  part  Rq,  the  frozen  fringe 
and  the  ice  layer  R2>  as  shown  in  Figure  1  whavrtj 
is  the  interface  between  Rj  and  Ri,  Hq  is  the  (PC 
isotherm  and  n  is  the  reference  boundary.  We  will 
assiune  that  file  pressure  of  water  is  kept  constant, 
Pn  at  n.  The  physical  properties  of  parts  Rq,  and  R2 
are  well  understood  but  our  knowledge  of  file 


Figure  1.  Schematic  drawing  of  a  steadUy 
growing  ice  layer  in  a  freezing^. 


physical  pn^rties  and  file  dynamic  behavior  of 
part  Ri  does  not  appear  sufficient  for  engineering 
applicatimis. 

Konrad  and  Morgenstem  (1980, 1981)  empiri¬ 
cally  found  that  fiie  mass  flux  of  water  /  at  file 
formation  of  the  so-called  final  ke  lens  is  propor¬ 
tional  to  the  average  temperature  gradient  in 
R}.  This  may  be  written  as 

SP^-f/Tf  (1) 

where  a  prime  denotes  dififerentiation  wifii  re¬ 
spect  to  x.  The  positive  proportionality  frulor  5P  is 
termed  the  segregation  potential,  ^lidi  is  a  prop¬ 
erty  of  a  given  s(^. 

Koirad  and  Morgenstem  (1982b)  also  found 
empirically  fiut  SP  depends  on  the  applied  pres¬ 
sure  Pa  and  fiw  pressure  of  water  Pq  at  iio  and  fiiat 
SP  is  a  decreasing  function  of  P,  and  file  sucticm 
-Po-  This  may  be  written  as 


SP=y(Po,PJ 

(2) 

—  y>0  ,  — y<  0. 
dPo  ap. 

(3) 

Extending  the  oono^  of  segregation  potential, 
Konrad  and  Morgenstem  (1982a)  fritroduoed  a 
semieir^iirical  mo^  (5P  nuxkl)  soil  freezing. 

Seve^  researchers  have  used  file  SP  model  to 
analyze  frost  heave  data  from  field  and  laboratory 
tests  (hfixon  1982,  IfriutssonetaL  1965,  Jessberger 
and  Jagow  1969).  Their  results  suggest  fiuit  the 
accuracy  of  file  SP  model  suffices  engineerit^ 
needs.  However,  the  SP  model  is  not  immune  to 
criticirtn.  Some  limitations  or  shortcomings  of  file 
SP  model  have  appeared  (bhizaki  and  Nishio 
1985,  Van  (Sassen  and  Sego  1989,  hfixon  1991). 


Nevertt\eles8,  it  is  important  to  point  out  that  an 
acxiirate  model  of  soil  freezing,  if  it  exists,  should 
be  able  to  explain  important  empirical  hndings 
such  as  eq  1, 2  and  3. 

Recently  the  steady  growth  condition  of  an  ice 
layer  was  studied  ma^ematically  and  experimen¬ 
tally  imder  a  negligible  applied  pressure  (Nakano 
1990,  Takeda  andNakano  1990,  Nakano  and  Take- 
da  1991)  and  under  an  applied  pressure  (Takeda 
andNakano  1993,  Nakano  and  Takeda  1994).  They 
have  shown  that  a  model  called  Mi  accurately 
describes  the  experimentally  determined  condi¬ 
tion  of  steady  growth. 

For  the  problem  of  a  steadily  growing  ice  layer, 
the  Ml  model  is  defined  as  a  ^zen  fringe  where 
ice  may  exist  but  does  not  grow,  and  where  the 
mass  flux  of  water /is  given  as 


f=-Kt^-K2^ 

for  z  in  Ri 

(4a) 

dx  dx 

K2/K1  -» y  as 

f-»o 

(4b) 

lim  P(3^  =  Pa 

X  -♦Ml 

(4c) 

X  in  Ri 

where  yis  a  constant,  P  is  the  pressure  of  unfrozen 
water  and  Kj  (i  =  1, 2)  is  the  transport  property  of 
a  given  soil  that  generally  depends  on  the  temper¬ 
ature  and  the  composition  of  the  soil.  The  Mi 
model  is  a  generaliution  of  somewhat  simpler 
models  (Derjaguin  and  Churaev  1978,  Ratkje  et  al. 
1982,  Kuroda  1985,  Horiguchi  1987)  in  which  die 
ratio  K2/K1  is  equal  to  y  regardless  of /. 

Experimental  methods  were  propo^  to  deter¬ 
mine  Ki  (Williams  and  Burt  1974,  Horiguchi  and 
Miller  1983)  and  K2  (Perfect  and  Williams  1980). 
According  to  Horiguchi  and  Miller  (1983)  Ki  of 
several  frozen  porous  media  is  described  in  the 
general  form  given  as 

Ki  =  AolTf*  (5) 

where  Aq  and  B  are  positive  material  constants. 
Nakano  and  Takeda  (1994)  experimentally  deter¬ 
mined  K2  of  Kanto  loam  in  the  form  given  as 

K2  =  /^“  A5T<0 

\#C2o|A/r|*  A>T 

where  K20  and  B  are  positive  constants  and  A  is  a 
small  negative  constant.  Since  empirically  deter¬ 
mined  functimts  and  X2  are  known  to  be  bound¬ 


ed,  the  functi(mal  form  of  eq  6  is  a  better  approxi¬ 
mation  to  their  actual  behaviors  in  a  neighborhood 
ofT  =  0“C. 

The  SP  model  is  semiempirical  in  nature.  An 
accurate  mathematical  model  is  needed  that  pro¬ 
vides  the  functional  dependence  of  SP  on  perti¬ 
nent  variables  specifying  given  thermal  and  hy¬ 
draulic  conditions  in  terms  of  well-defined  func¬ 
tions  (or  parameters)  describing  the  properties  of 
a  given  soil.  In  this  report  we  will  present  the 
results  of  our  study  on  the  problem  of  a  steadily 
growing  ice  layer  by  using  the  Mj  model.  We  will 
show  that  the  Mj  model  provides  the  functional 
dependence  of  SP  that  is  coiuistent  with  empirical 
findings  (eq  1, 2  and  3). 

PROPERTIES  OF  Ml 

An  analytical  solution  was  derived  to  the  prob¬ 
lem  of  a  steadily  growing  ice  layer  imder  the  Mi 
model  (Nakano  1990).  A  schematic  drawing  of  a 
typical  temperature  profile  is  shown  in  Figure  2, 
where  the  temperature  is  continuous  but  the  tem¬ 
perature  gradient  is  discontinuous  at  rii  where  the 
ice  layer  grows.  Since  the  amount  of  heat  trans¬ 
ported  by  convection  is  generally  much  less  than 
that  by  conduction,  the  temperature  profiles  in  R^, 
and  R2  in  the  vicinity  of  Ri  are  nearly  linear  (Na¬ 
kano  and  Takeda  1991). 

The  temperature  T]  and  the  limiting  value  of 
the  temperature  gradient  at  tti  are  given  approxi¬ 
mately  (Nakano  1990)  as 


Ti  =  — a  ao 5 

(7a) 

T'(n+)  =  -boo 

(7b) 

where  Oq  is  the  absolute  value  of  the  temperature 
gradient  at  Hq,  5  is  the  thickness  of  Ri  and  de¬ 
notes  the  limiting  value  as  x  approaches  while  x 
is  in  R].  The  a  and  b  in  eq  7a  and  7b  are  correction 
factors  that  take  account  of  effects  of  the  convec¬ 
tive  heat  transport  and  the  variable  thermal  con¬ 
ductivity  in  Ri  (Nakano  and  Takeda  1991).  When 
these  effects  are  negligible,  a  and  h  are  equal  to  one . 
We  will  neglect  these  ejects  hereafter. 

When  an  ice  layer  is  steadily  growing,  accord¬ 
ing  to  the  Ml  model,  ice  does  not  grow  in  R^;  the 
growth  of  ice  occurs  only  at  nj.  The  balance  of  heat 
for  R]  is  given  as 

ktai  =  kooo  +  (L  -  C2  Ti)/  (8) 

where  Jt]  and  Ag  are  the  thermal  conductivities  of  R2 


and  L  is  the  latent  heat  of  fusion  of  water,  and 
C2  is  the  heat  capacity  of  R2.  is  the  absolute  value 
of  dte  limiting  temperature  gradient  at  defined 
as 


the  formatiiHiof  the  final  ice  lens.  For  given  a  and 
Sq,  the  necessary  and  sufficient  condition  for  the 
steady  growfii  of  an  ke  layer  is  given  (Nakano  and 
Takeda  1991,  Nakano  and  Takeda  1994)  as 


ai  =  -  lim  t'W  =  (9) 

xinRi 

wherenfdenotes  delimiting  value  asxapproadies 
n-i  while  x  is  in  R2- 

A  schenutic  drawing  of  a  typical  pressiue  pro¬ 
file  is  given  also  in  Figure  2  where  P„  is  kept  at  0.1 
MPa.  The  pressure  of  water  P  is  continuous  in  the 
combined  region  of  Rq,  and  Rj  but  the  gradient  of 
P  may  be  disomtinuous  at  hq.  Darcy's  law  clearly 
holds  true  in  Rq.  Neglecting  ttie  gravity  efiect,  we 
obtain 

AP  =  P„-Po=WV  (10) 

where  =  no  -  n  >  0. 

Ithas  been  shown  (Nakano 1992)  thataccording 
to  the  Ml  model,  the  segregation  potential  func¬ 
tion  y  is  given  as 

y  =  -/vr(m+)  =r  =  KzM  (11a) 


(Jti/ito)oci>ao2iti(ito-«'Lyr'cci-  (Hi’) 

The  condition  of  eqllbis  explained  in  de  diagram 
oftemperature  gradients  (Rg.  3).  inRgureSis 

given  as 

o(o  =  (ki/jb))ai  (m  RT*.  (11c) 

If  an  existing  ice  layer  neither  grows  nor  melts,  / 
vanishes  cm  Iff  *.  We  will  refer  to  the  regirm  as 
where  Oq  >  <  0  arui  an  ice  layer  is  melt¬ 

ing.  iff  is  given  as 

00  =  lfci(lfco  + Ly)"*ai  cm  Rt  .  (lid) 

The  bcnmdary  iff  is  a  curve  sterruning  from  de 
origin  becauseydepends  cm  Oq.  The  steady  growd 
of  an  ice  layer  occurs  in  de  regicm  P,.  The  forma¬ 
tion  of  de  final  ice  lens  occurs  on  Iff  where  eq  11a 
holds  true. 

The  temperature  Ti  is  given  (Nakano  1990, 
Nakano  and  Takeda  1^)  as 


3 


in 


Figure  3.  Diagram  of  temperature  gradiertts  Oj  and  Oq. 
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Figure  4.  Sdtematic  drawing  of  a  steadily 
growing  ice  layer  under  an  applied  load. 


-Ti  =  (o  +  W^// 

(12a) 

/  *  Y^ai-RoV®o* 

(12b) 

where  0,  and  are  defined  as 

o  =  P.-P„ 

(13a) 

♦01  (Ti)  =  77*  I  (Xo/Ki)  ('X2/X2o)dT  (13b) 

♦m(Ti)  =  77*J  (Ko/Ki)dT. 

(13c) 

Using  eq  1  la,  12a  and  12b,  the  mass  flux/*  at  the 
formation  of  the  final  ice  lens  is  given  as 


f*=(^r<y8o)[-rf/-oi  (i4a) 

linUy^axM-Ko'^KiinhuiTf).  (14b) 

Equatimis  14a  and  14h  imply  that/*  is  uniquely 
determined  by  T*  when  and  o  ate  given.  Using 
eq  11a,  we  will  r^uce  eq  14a  and  14b  to: 

K2m  =  Aro(oo8o)-«  [-7T  /-ol.  (15) 

Since  X2  is  an  increasing  function  of  T f ,  and  Tf 
are  one-to-one  and  Tf  may  be  ccmndeted  to  be  a 
functiai  of  K2.  Therefore,  eq  15  implies  thaty  is  a 
function  of  Og,  Sg  and  o. 


y  =  y(okv8ivO)-  (16) 

Ithasbeen  shown  empirically  (RaddandOertle 
1973,  Takashi  et  al.  1981)  that  there  is  a  unique 
temperature  T**  at  for  a  given  o  >  0  when  an 
existing  ice  layer  neifoer  grows  nor  melts  and  / 
vanishes.  This  temperature  Tf*  at  the  phase  equi¬ 
librium  of  water  is  given  as 

a^-yrr,  if  /=0.  (17) 

When  an  ice  layer  is  steadily  growing  imder  a 
positive  a,  is  less  than  T**  (Nakano  and  Takeda 
1994)  and  the  frozen  fringe  Ri  may  be  divided  into 
two  regions  Rio(0  ^  T ^  r?*)  and  Rii(rr*  >  T  > 
ro  as  shown  in  Figure  4.  At  the  phase  equilibrium 
/  vanishes  and  Rn  also  vanishes.  From  eq  4b  we 
obtain 

K2/Ki=Y.  in  Rio  if/=0.  (18) 

It  is  easy  to  see  that  eq  12a  is  reduced  to  eq  17  at  the 
phase  equilibrium  because  of  eq  18.  An  important 
question  arises  wheflier  or  not  eq  18  holds  true 
when  /  does  not  vanish.  Since  fQ  (t  =  1, 2)  mainly 
depends  on  T  and  the  composition  of  R^q,  it  is 
prbbable  that  eq  18  holds  true  for  /  >  0  unless  / 
significantly  aff^  the  composition.  We  will  as¬ 
sume  the  v^dity  of  eq  18  regardless  off.  With  this 
assumpticm  eq  4a  beeves  one  step  closer  to  flte 
simpler  models  described  earlier. 
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FUNCTION  y 

Xg  =  Ki/b2 

(24e) 

The  functional  form  of  y  depends  on  those  of  Xj 
and  X2.  Based  on  available  experimental  findings 

Xj  =  bj — 1>2  1* 

(24f) 

(Horigudii  and  Miller  1983,  Nakano  and  Takeda 
1994),  we  will  assume  that  Kj  and  X2  are  given  as 

For  a  special  case  where  X| 

= 0,/]  and/are  given  as 

Xi  =  (^®  AST<0 

/i  =  -l^ln(y/X2o) 

(25a) 

lXo(A/T)^  A>T 

/(0)  =  -fciI^ 

(25b) 

X2  =  /^“  AST<0 

A<  Tf* 

(25c) 

\K2o(A/tP  A>T 

\-bi*  +  ln2 

A^Tf* 

Because  of  eq  4b  Kq  and  Kjo  are  related  as 

X2o/Ko  =  Y.  (21) 

The  values  of  parameters  in  eq  19  and  20  deter¬ 
mined  empirically  (Nakano  and  Takeda  1994)  for 
Xante  loam  are 

Xo  =  1.77  X 103  g/(cm  d  MPa) 

Xjo  =  1.98  X 103  j  o^) 

=  0520 
b2=1.04 

v4  =  -1.5xl0-*“C 


where  d  denotes  day. 

Using  eq  19, 20  and  21  and  eliminating  Tt 
reduce  eq  15  to 

*,we 

/i  -  /(o)  =  «o  So  Co  (1  +  /2)y 

(22) 

where 

/i  =  Ciy-*<' 

(23a) 

c,=X7'x^Mfri+i) 

(23b) 

7(0)  =  x;'  (62- *1) 

(23c) 

Xi^  (62  —  bi)  +  z  A  <  Ti 

/(a)  = 

Xi  2  A^T\ 

(23d) 

Co  =  -[AKmP>0 

(23e) 

/2  =  C2(ao8o)‘^  >  0 

(24a) 

C2  =  - Abi(bi  +  l)-^>0 

(24b) 

2  =  C30 

(24c) 

C3  =  -(YAr'>0 

(24d) 

It  is  clear  from  eq  10  and  11  diat  Oo  ^‘‘^d  ^  must 
be  positive  because  die  mass  flux/*  is  positive  at 
die  formadm  of  die  final  ice  lens.  Whoi  OgSg  is 
very  small,  eq  22  is  reduced  to 

/i-/{o)  =  CoC2y  (26) 

It  follows  from  eq  26  diat  y  no  longer  dqiends  on 
Oo  or  Sq  when  is  voy  smalL  As  shown  in 
Appendix  A,  eq  22  generally  possesses  a  unique 
positive  rooty.  However,  under  certain  conditions 
this  unique  root  may  not  be  physically  meaning¬ 
ful.  This  implies  that  die  formation  of  die  final  ice 
lens  may  not  take  place  under  sudi  conditims. 

When  eq  22  possesses  as  a  meaningful  root,  we 
will  study  ^e  dependence  of  y  on  parameters,  Og, 
§0,  AP  and  obelow.  We  will  write  eq  22  as 


Ji  -  J{o)  =  ao5oCoy  +  Co  Ciy.  (27) 

Diflerentiating  eq  27  with  respect  to  (Xq  and  fig,  we 
obtain 

f-5oCoyli~*  A,i^0  ^28a) 

\-6oCoyUo'  ^1  =  0 

/-otoCoyU”^  ^28b) 

^  \-aoCoylio'‘  Xi  =  0 
where 

U  =  otoSoCo  +  C0C2  +  CiXoy-(^*»)  (29a) 

Uo  =  ao8oCo  +  C0C2  +  b2  (29b) 

CiXo=X^/(ftx  +  l)>0.  (29c) 


It  follows  from  eq  28a  and  2%  that  y  is  a  decreasing 
functicai  of  both  ocg  and  5g. 
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From  eq  10  and  11,  we  obtain 

ooSoCoy  ^  3.37oco6oy 

(35e) 

aoSo=  KoAf^“*. 

(30) 

CoC2y  =  1.75xlO~‘y 

(35f) 

Using  eq  30,  we  write  eq  22  as 

Ji  -  J(o)  =  KoCoAP  +  CoC2y.  (31) 

Differentiating  eq  31  with  respect  to  AP,  we  obtain 

I-KoCoKT'  X,i  =  0 

where 


where  units  of  variables  are  y  [g(cm  'C  d)"*],  a 
[MPa],  Oq  I*C  cm"*]  and  6o  [cm]. 

It  is  easy  to  find  from  eq  35e  and  35f  that  die 
term  CgCjy  is  much  less  than  the  term  OgSoCoy 
unless  OgSo  is  very  small  Neglecting  this  small 
term,  we  reduce  eq  27  to 

h  -  Kp)  -  ctoSoCoy.  (36a) 

Using  AP,  we  write  eq  36a  as 


V»CoC2+CjXoy-<*«'*'l  (33a) 

Vb=CoC2  +  6r'y-'.  (33b) 

It  follows  from  eq  32  that  y  is  a  decreasing  function 
of  AP.  This  property  of  y  is  consistent  with  the 
empirical  finding  (the  first  inequality  of  eq  3). 

Differentiating  eq  27  with  respect  to  o,  we  ob¬ 
tain  for  X]  ^0: 


/i-/(o)  =  KoCoAP  (36b) 

When  o  =  0,  the  term  /(O)  is  much  less  than  the 
other  two  terms  of  eq  3te  or  36b.  When  /(O)  is 
neglected,  eq  36a  and  36b  are  reduced  to 

y  =  (Cl  (ao5o)‘‘'^  *  (37a) 

y  =  [Cl  /  {KoCo)f^  (A P)"’'^  (37b) 


^^(-C3U-‘ 

\-C3^V*‘^U'‘ 

When  Xi  =  0,  we  obtain 

f^^/-C3Uo' 


A  <  Tf* 
A  S  Tf* 


.(34a) 


A  <  Tf ' 
A  S  TV 


.(34b) 


It  is  easy  to  find  from  eq  34a  and  34b  that  y  is  a 
decreasing  function  of  a.  This  property  of  y  is  con¬ 
sistent  wi  A  the  empirical  finding  (the  second  ine¬ 
quality  of  eq  3). 

Using  experimentally  determined  parameters 
of  Kanto  loam  (Nakano  and  Takeda  1994),  we  will 
show  the  behavior  of  y  of  Kanto  loam  below.  For 
Kanto  loam  the  values  of  ^  and  Xi  are  0.462  and 
0.480,  respectively.  We  will  present  the  actual  val- 


ues  of  ea^  term  in  eq  27: 

/i  =  4.75xl0y-®-^ 

(35a) 

7(0)  =  1.08 

(35b) 

2  =  5.95x103® 

(35c) 

7(0)  = 

|l.08  +  5.95x10*0 
16.48  X  10  o®-^ 

^<^‘(35d) 

A^TC 

Suppose  that  we  measure  a  set  of  values  either  (y, 
Oj,^)  or  (y,  AP)  and  that  my  and  In  (or  In  AP) 

are  linear.  If  we  know  one  of  three  parameters,  bi, 
b2  and  A,  the  remaining  two  parameters  can  be  de¬ 
termined  by  eq  37a  (or  eq  37b).  When  ct  ^  0,  the 
term  /(o)  is  not  generally  negligible  and  we  must 
use  eq  36a  (or  eq  36b)  to  determine  an  accurate 
dependence  of  y  on  OoSq  (or  AP). 

In  Figure  5  we  plotted  experimental  data,  y  vs. 
Oq  with  = 2  cm  imder  o = 0, 16.2, 48.7  and  195  kPa 
(Nakano  and  Takeda  1994)  together  with  predict¬ 
ed  y  calculated  by  eq  36a.  It  is  easy  to  see  that  In  y 
and  In  Oq  are  nearly  linear  when  o  =  0.  As  a  in¬ 
creases  the  relationship  between  In  y  and  In  Oq 
becomes  nonlinear  as  anticipated. 

CONCLUDING  REMARKS 

Many  models  of  ice  segregation  have  been  pro¬ 
posed  in  the  past.  However,  the  SP  model  pro¬ 
posed  by  Konrad  and  Morgenstem  (1980, 1981)  is 
(me  of  few  that  were  built  cm  an  empirical  base. 
The  SP  model  has  been  ccmveniently  used  to  solve 
engineering  problems.  However,  we  need  an  ac¬ 
curate  mathematicalmodel  ffiatprovides  the  func¬ 
tional  dependence  of  SP  (m  pertinent  variables 
specifying  given  thermal  and  hydraulic  condi¬ 
tions  in  terms  of  well-defined  functions  (or  para¬ 
meters)  describing  the  properties  of  a  given  soil. 
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Li  response  to  such  a  need  Nakano  (1990)  pro¬ 
posed  a  new  madiematical  model  called  M|.  Ef¬ 
forts  have  been  made  to  validate  ttie  model  by 
empirical  findings  and  experimental  data  (Takeda 
andNakano 1990,  Nakano  andTakeda  1991,  Take¬ 
da  andNakano  1^3,  Nakano  and  Takeda 1994).  In 
this  report  we  presented  foe  result  of  our  efforts  to 
validate  foe  model  by  using  empirical  findings 
foat  were  used  to  build  the  SP  model.  It  is  impor¬ 
tant  to  mention  that  eq  1, 2  and  3  are  empirical  find¬ 
ings,  not  assumptions.  If  foe  Mj  model  is  accurate, 
it  must  be  able  to  explain  these  important  em¬ 
pirical  relationships.  We  have  shown  that  foe  Mi 
model  can  explain  these  relationships. 
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APPENDIX  A:  POSITIVE  ROOT  OF  EQUATION  22 

We  will  show  that  eq  22  possesses  a  unique  positive  root  y  >  0.  Using  eq  27  that  is  equivalent  to  eq 
22,  we  will  introduce  a  fur^on  F(y)  dehn^  as: 

F(y)  =  /i  -  «o6oCo  y-  C0C2  y  -  /  (a).  (Al) 

It  is  clear  that  a  root  of  eq  22  satifies 

F(y)  =  0  (A2) 

Differentiating  F(y)  with  respect  to  y,  we  obtain 

F'(y)  =  J[(y)  -  (aoSoCo  +  C0C2)  (A3) 

where  a  prime  denotes  the  differentiation  with  respect  toy. 

From  eq  23a  and  25a  we  obtain 

^'*0  (A4) 

Since  and  ^2  are  positive,  /'(y)  is  negative.  Hente  r(y)  is  also  negative,  in  other  words,  F(y)  is  a 

strictly  decrea^g  function  of  y.  We  will  examine  the  b^vior  of  F(y)  for  two  cases:  Case  1,  Xj  2: 0, 
and  Case  2,  Xj  <  0. 

For  Case  1,  F(y)  approaches  +«» as  y  approaches  0  while  F{y)  approadies  -<>o  as  y  approaches  +». 
Since  F(y)  is  a  strictly  decreasing  function,  eq  22  possesses  a  uiuque  positive  root.  For  Case  2,  F(y) 
approaches -/(o)  as  y  approaches  0  while  F(y)  approaches -o®  as  y  approaches  +«».  WhenXj  <0,b2> 
bi  +  1.  Hence  /(a)  is  negative  and  eq  22  possesses  a  unique  positive  root.  It  may  be  noted  paren¬ 
thetically  that  A  <  Tf*  implies  2  <  1  because  of  eq  17 

We  have  shown  that  eq  22  possesses  a  unique  positive  root  regardless  of  Xj.  When  an  ice  ]ayes  is 
growing,  the  temperature  must  be  less  than  Tf  *.  Ihis  implies  that  y  must  be  less  than  Ki 
namely: 

y  <  y*  *  =  fC2(ri  •)  =  K20 (AS) 

It  should  be  noted  that  y**  =  K20  when  o  =  0.  Equatiwi  A5  implies  flat  F(y**)  must  be  negative;  that 
is 


F(y**)  =  /3(y*’,a)-C42-*»<o 

(A6) 

where 

/3  =  /<y**)-/(o) 

(AT) 

C4  =  X2oCo(o(o8o  +  C2)  >  0. 

(A8) 

Ihe  functicm  is  given  as: 
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1bibi'>0, 
lm*b\bi^-z 
bl'  >  0, 


Xi  ■  0 ,  o  ■  0 
Xi»0,A<7f  • 

X\^0,Akti* 


(A9) 


[fri  bi  >  0, 


Xi  =  0 ,  o  =  0 


/3  = 


Inz  +  61^2^ 

]b2^  >  0, 


z 


Xi  =  0,A<  If 
Xi  =  0 ,  A  2  Ii* 


(AlO) 


It  follows  from  eq  A6,  A9  and  AlO  that  the  condition  of  eq  A6  is  always  satisfied  when  Xx  ^  0  and  A 
2  n* .  However,  for  oliter  cases  it  is  not  certain  ttiat  eq  A6  is  satisfied. 
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